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Abstract 

We consider the problem of multiparticle production at threshold in a (/)^-theory 
with an 0(A'^i+A''2) symmetry softly broken down to 0{Ni) x 0{N2) by nonequal 
masses. We derive the set of recurrence relations between the multiparticle ampli- 
tudes which sums all relevant diagrams with arbitrary number of loops in the large- 
A'^ limit with fixed number of produced particles. We transform it into a quantum 
mechanical problem and show how it can be obtained directly from the operator 
equations of motion by applying the factorization at large A^. We find the exact 
solutions to the problem by using the Gelfand-Dikii representation of the diagonal 
resolvent of the Schrodinger operator. The result coincides with the tree amplitudes 
while the effect of loops is the renormalization of the coupling constant and masses. 
The form of the solution is due to the fact that the exact amplitude of the process 
2^n vanishes at n>2 on mass shell when averaged over the 0(A^i^2)-iiidices of in- 
coming particles. We discuss what dynamical symmetry is behind this property. We 
also give an exact solution in the large- iV limit for the model of the 0{N)+ singlet 
scalar particle with the spontaneous breaking of a reflection symmetry. 
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1 Introduction 

The problem of calculating amplitudes of multiparticle production at threshold has re- 
cently received a considerable interest [jl|-||2^. The explicit results |^, for the tree level 
amplitudes in the 0^ theory demonstrate the factorial growth which is expected due to the 
large number of identical bosons in the final state. While the threshold amplitudes have 
a vanishing phase volume, they are a lower bound for tree amplitudes with nonvanishing 
spatial momenta of produced particles which violate unitarity [^ H, ^. Therefore, loops 
should be taken into account. 



The one-loop calculation of the l^n amplitude discovered [|T2| a very interesting 



property of the tree 2— i>n threshold amplitude which is that it vanishes for n>4 when 



incoming particles are on mass-shell. This nullification has been extended ||T3[-p5[ to 
more general models and is crucial for calculations of the amplitude 1— i>n at the one- 
loop level. A dynamical symmetry which may be responsible for the nullification has 



been discussed in Ref. [Q. An interesting question is which properties of the tree and 
one-loop amplitudes could survive in the full theory which includes all loop diagrams. 

The simplest exactly solvable model is the 0(iV)-symmetric 0"^ theory in the large- A^ 
limit which is solvable |^ since only the bubble diagrams contribute if there is no multi- 
particle production. In Ref. |2^ the set of recurrence relations between the amplitudes of 
multiparticle production at threshold has been derived in the large- A^ limit at fixed num- 
ber of produced particles n which sums all relevant diagrams with an arbitrary number 
of loops. These recurrence relations have been transformed into a quantum mechanical 
problem whose exact solution has been found through the use of the Gelfand-Dikii rep- 
resentation of the diagonal resolvent of the Schrodinger operator. The result was quite 
similar to the tree amplitudes while the effect of loops is the renormalization of the cou- 
pling constant and mass. The form of the solution was due to the fact that the exact 
amplitude of the process 2— >n vanishes for n>2 at large A^ on mass shell when aver- 
aged over the 0(A^)-indices of incoming particles. This nullification occurs for dynamical 
reasons because of the cancellation between diagrams. 



In the present paper we extend the approach of Ref. |£6| to more complicated models. 
We consider the model of two scalar fields: a A^i-component field with the mass mi and 
a A^2-component one with the mass m2- The interaction is 0(A^i+A^2) invariant while 
the mass term breaks this symmetry down to 0{Ni) x 0{N2) for nonequal masses irti 
and 1712- We consider also the case when the second field is a single component one 
{i.e. N2 = 1) and has a wrong sign of the mass square, tti^ < 0, which results in the 
spontaneous breaking of a reflection symmetry. The motivations for studying this model 
are as follows: 

1. The tree amplitudes 1—^n at threshold can be found for the model with the softly 



broken 0(A^)-symmetry extending the results of Ref. pSf given for N=2. A question 
arises as to how these results, obtained exploiting in analogy with 1-dimensional 
integrable systems, can survive at large A^ in higher loops. 



2. The model with A^2=l and m| < describes after the spontaneous breaking of 
a reflection symmetry 1 particle with mass A/2|m2| and A^— 1 particles with mass 



ml + |r?7,2|. If the 0(A^)-symmetry were not explicitly broken by masses {i.e. if 
ml = ml), these A^— 1 particles would be massless — the Goldstone bosons as- 
sociated with a spontaneous breaking of the 0(A^)-symmetry — so that the con- 
sideration would not be possible. Having the large-A^ solution, one can perform 
a comparison with the one for A^=l which was obtained by semiclassical meth- 



ods 16, 22 . 



3. The model with the softly broken 0(2)-symmetry reveals at the tree level a dy- 
namical symmetry which is responsible for the nullification in the case of further 
restrictions on kinematics when the mass ratio is an integer so that no spatial mo- 



mentum is involved in the decay 1-^n [0]. The problem is whether this dynamical 



symmetry exists at arbitrary rrii and m2 even at the tree level. 
We derive the set of recurrence relations between the multiparticle amplitudes in the 
large- iV limit at fixed number of produced particles. We transform it into a quantum 
mechanical problem and show how it can be obtained directly from the operator equations 
of motion by applying the factorization at large A^. We find that the quantum mechanical 
problem can be exactly solved through the use of the Gelfand-Dikii techniques. Thus we 
get the generating functions for the amplitudes where the effect of loops is reduced to a 
renormalization of the coupling constants and masses. An interesting new phenomenon 
that shows up in the 0{Ni) x 0{N2) model is an appearance of an imaginary part in 
the amplitudes which is due to inelastic processes. An extension of that result to the 
case of 0{Ni) x ... x 0{Ns), s > 2, is shown to be straightforward. For the model of the 
0{N) + singlet scalar particles, the generating function for the amplitudes in the large 
A^ limit (no loops of a singlet field should be taken into account) can be extracted from 
the explicit solution for the 0{Ni) x 0{N2) model. 

An interesting modification of the result appears when we consider a model with a 
negative mass square of the singlet particle. This model corresponds to a spontaneous 
breaking of the parity reflection of the singlet field. The amplitudes turn out to be real for 
this model. By using the exact solution for this problem we demonstrate a nullification 
of amplitudes except for the case of production of 2 0{N) particles and 1 physical singlet 
particle. 

This paper is organized as follows. Sect. 2 is devoted to the definitions and the 
description of the kinematics. In Sect. 3 we derive the set of recurrence relations for the 
0{Ni) X 0(A^2) model and transform it into a quantum mechanical problem. In Sect. 4 we 
present an exact solution and demonstrate that it satisfies the set of equations. The exact 
solution is possible since the diagonal resolvent of the Schrodinger operator has a very 
simple form. In Sect. 5 we consider the 0{N) + singlet particle model with spontaneously 
broken parity reflection. We conclude by summarizing the results of this paper and give a 



generalization of the dynamical symmetry of Ref. [53] to the 0{Ni) x 0{N2) type models 



2 Definition of the amplitudes 

We consider a model of an 0{Ni) scalar field multiplet x" ? ct = 1, ■■■,Ni , coupled to a 
0{N2) scalar field 0", a = 1, ..., A^2- The theory is defined in the 4-diniensional Minkowski 
space by the following Lagrangian 

£ = lid.x'mx') + Ud,r){d,r) - imiixV) - imiirr) - iHrr + xVy (2.1) 

where the summation over repeated indices is implied. The interaction term is 0{Ni+N2) 
invariant while the mass term possesses the 0{Ni)xO{N2) symmetry for nonequal masses 
rui and m2. 

We denote by Al_^ f,^^^^ ^^{n, k) the amplitude of production of n on-mass-shell parti- 
cles X with the ©(A^^i) -indices bi . . . bn and of k on-mass-shell particles (p with the 0{N2)- 
indices ai, . . . , a^ at rest by a (virtual) particle x with the 0(iVi)-index b and the energy 
nrrii + km2'- 

K...b„,a,...aS^, k) = {{nrm + km2f - mDm^'^m^'' {b^ . . . bn, «i • • • «fc|x'(0)|o) . (2.2) 

We also denote by BJI^ ^^ ai...aX^i k) the amplitude of production of n on-mass-shell par- 
ticles X with the O(A^i) -indices 61 . . . 6„ and of k on-mass-shell particles with 0{N2)- 
indices ai, . . . , a^ at rest by a (virtual) particle with the 0(A^2)-index /5 and the energy 
nrrii + km2'- 

K...K,a,...a, in, k) = ((nmi + fcms)' - mDm^^m:^^ (b, . . . 6„, a,... «fc|0^(O) |o) . (2.3) 

The extra factors in Eqs. (|2.2| ) and ( p.3|) are due to the amputation of the external lines. 

It is convenient to multiply these amplitudes by the vectors C,i and ^2^ which results in 
a symmetrization over the 0{Ni) and 0(A^2)-indices of the produced particles and define 



An, k) = <...,„,«,..„, (n, k) e^' . ..e,"^r ■ --iT , (2.4) 

and 

5"(n, k) = 5,"„,.,„,,,...„,(n, k) e^' ■ --ei-^r ■ --^r ■ (2.5) 

The amplitudes A'^ln, k) and B°'{n, k) are depicted graphically in Fig. 0a, b. 

The other quantities which appear in the set of the recurrence relations between the 
amplitudes are the amplitudes Df^ j^^^^ ^ {n,k;p) and D'^^,^^ai...aXn,k]p) of the pro- 
cesses when two particles a and b (or a and j3) with the 4-momenta p + nqi + kq2 and —p, 
respectively, produce n on-mass-shell particles x with the 0(A''i)-indices 61,..., 6„ and 
the equal 4-momenta gi = (rrii, 0) and k on- mass-shell particles with the 0(A'"2)-indices 
«!,..., ctfc and the equal 4-momenta g2 = (^^2, 0) in the rest frame. We define again 

D?(n, k- p) = <..,„,«„..„, (n, k; p) i\' ■ ■ ■ il^ir " " " iT , (2-6) 



A^'in.k) = 




B'^in.k) = 

b 




a) 



b) 



Figure 1: The graphic representations of the multiparticle amplitudes A"(n, k) (a)) and 
B'^{n,k) (b)). 
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Figure 2: The graphic representations of the multiparticle amplitudes D'^{n,k;p) (a)) 
andZ?"^(n,A;;p) (b)). 



and 

Df{n, k- p) = <„,„,„,...„, in, k; p) t^' ■ ■ • ^t^!^^ ■ ■ ■ ^f • (2.7) 

The amplitudes D'^(n,k;p) and D" (n,k;p) are depicted graphically in Fig. |^a, b. The 
definitions (|2.4|), (p.5|), ( |2.6| ) and ( pTT] ) generalize those of Ref. [^] to the case of nonequal 
masses mi and m2- 

Notice that we do not introduce the amplitudes associated with the mixed case when 
the incoming particles are x ^^^ 0- This is because the unsuppressed in the large N limit 
loop diagrams contain only loops of Xb and 0q, particles with contracted 0(iVi 2) indices 
(it is assumed that A ~ 1/^)- Notice also that in the case Ni = N ^ 00 and N2 =finite 
we can simply drop out the contributions from the loops with virtual particles. That 
means that in this case we have to consider only D"^^ vertices. Actually we shall see that 
the amplitudes for the Ni = N ^ 00, N2 =finite case can be easily extracted from 
the expression for the amplitudes given below for the case A^i,2 ^ 00 hj taking formally 
N2/Ni = 0. 

It is easy to estimate at the tree level that 

A\n,k) ~ A^^ (e?)^ {a^'^e,, B^in,k) ~ A^^ (e?)^ {^i)"^ ^2, (2.8) 
where ^f and ^2 stand for ^"^" and ^f ^^ respectively. Since 

A ~ 1 (2.9) 



in the large- iV limit p8|, we choose 

^^2 _ ^2^ _ A^ (2.10) 

for all the amplitudes ( p.8| ) to be of the same order in 1/A^. 



3 Schwinger— Dyson equations at large N 



3.1 The diagrammatic derivation 



The condition (|2.9| ) leaves at large- A^ only the diagrams with the largest possible number 
of sums over internal ©(A"! 2)-indices which propagate along closed loops of diagrams. 
The proper recurrence relations which extend the usual Schwinger-Dyson equations to 
the case when n particles x and k particles are produced are depicted for A^{n, k) and 
B°'{n,k) in Fig. |^. The solid lines are associated with the propagation of the 0{Ni)- 
indices (the field x), while the dashed ones stand for the propagation of the 0(A'2) indices 
(the field (p). Each vertex is the sum of three possible permutations of the 0(A^)-indices. 
This is taken into account in Fig. ^ by a combinatorial factor. The graphic notations 
becomes clear if one introduces the auxiliary field a{x) = x'^i^) + (p'^{x) which propagates 
in the empty space inside vertices. 
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+ 




Figure 3: The recurrence relation for the amplitude A'^(n,A;). The analytical formula 
is given by Eq. ( |3.1| ). The analogous recurrence relation for the amplitude 
B"{n, k) can be obtained by interchanging the solid and dashed lines. 



The recurrence relations of Fig. ^ read analytically as follows 

n\k\ 

A [n, k) = A 2_^ On.ni+n2+n3^k,ki+k2+k3] 



' ni\n2\n3\ki\k2\ks\ 
A-{n,,k,)A\n2,k2)A\ns,ks) 






{{nirrii + ki'm2Y ~ m\){{n2mi + k2m2Y ~ m\){{n^ini + A;3m2)^ — ml 

+ ^ / , (^n^ni+n2+n-j,0k,ki+k2+k'i 



+ 



ni\n2\n'j,\ki\k2\kz\ 
A%n,,k,)B''{n2,k2)B"{n^,ks) 



n-^,k2 .ko^odd. 
no ,nq ,fci ^even 






+ 



ni!n2!A;i!A;2! 

n2 ,fc]^ ,^2— ei'eri 



/7^KV,A:2;p) + /^r(^2,fc2;p)) , (3.1) 



((riiTTii + kim2Y ~ ml) J (27r)^ 
and similarly 

n\k\ 



B [n, k) — A 2_^ On,ni+n2+n30k^ki+k2+k3- 



fei.„2.n3=odd, • ni\n2Wki\k2\ks\ 



n^ ,^2 ,fc3^ei;en 

B''{n,,k,)A\n2,k2)A\n3,k:, 



{{nirrii + kim2Y — m2){{n2m,i + k2m,2)^ — mf)((n3mi + k^m2Y ~ ^1 



+ 



ni!n2!n3!A;i!A;2!A;3! 
B^{n^,k^)B^{n2,k2)B^{n^,k:,) 



k^^odd. 



{{niiTLi + kim2Y ~ Tn'2){{n2Tni + ^2^2)^ ~ '^Dll^s^^i + ^3^^2)^ ~ ^2 

+ -^ Z^ 0n,ni+n2'^fc,A;i+fc2 



+ 



ni\n2\ki\k2\ 



ni .no , fco— euen 



/7^KV2,fc2;p) + /^rKfc2;p)) . (3.2) 



((niTTii + kim2Y ~ m,2) J (27r) 

Notice that only the summed over the 0{Ni) and 0(A^2)-indices quantities D^{n2, k2]p) 
and D2°'{n2, k2]p) enter these recurrence relations. 

It is the property of large A'^1^2 that the recurrence relations for D'^{n,k\p) and 
D'^{n,k]p) express it via A^, i?" and D^, D^" again. The recurrence relations for 
D°^°-{n, k]p) and D^^iji^ k;p) are depicted in Fig. H. 



n, k 




+ 





\ri2,k2) 
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+ 




Figure 4: The recurrence relation for the amplitude D'^(n, k;p). The analytical formula 
is given by Eq. ( |3.3D . The recurrence relation for the amplitude D2'^{n,k;p) 
can be obtained by interchanging the solid and dashed lines. 



The recurrence relation of Fig. ^ reads analytically as follows^ 






1 A\ni,k{)A\n2,k2)Dl\n^,k^;p) 

+ 

ij) + nagi + ^352)^ - "^1 ((^i"^i + kim2)'^ - m\){{n2mi + k2m2y - m\) 

1 5/3 (m, fci)S/5(n2, k2)Dl\n^, k^; p) 

~l~ 

{p + n^qi + ^3^2)^ - "^1 (("-i^^i + kim2Y - ml){{n2mi + A;2"^2)^ - "^2) 

+ A 2^ 5n,ni+n2^fc,A:i+A:2 j ,, ,, , 

ni,ki=even x i x ^ 

1 /• d^A; 



/ 72^^^''^'''' ^'' ^^ + ^^^^^""'^ ^'' ^))^r(^2, ^2; p) . (3.3) 



(p + n2qi + ^252)^ - ml J (27r)'' ' 
Similarly we have 

-^0 {n,k)=X 2_^ On,ni+n2+n30k,ki+k2+k3- 



n„n2=odd, ■ ni\n2Wki\k2\ks\ 

no ,k^^even 

1 A\n,,k,)A\n2,k2)D^''{ns,ks;p) ^ 

{p + n3gi + ^352)^ - ^2 {{nimi + kim2Y - mi)((n2mi + k2m2Y - mf) 

+ A 2^ ^n,ni+n2+n3^fc,A:i+fc2+fc3 , , ,, ,, ,, , 
fci,fc2=odd, ni!n2!n3!fci!fc2!/C3! 

1 B''{n^,k^)B%n2,k2)D%''{n^,k^-p) ^ 

{p + n^qi + k-^q2Y - ml {{uimi + kim2Y - m2){{n2mi + k2m2Y - m^) 

~\~A 2_j ('■n,ni+n2^k,ki+k2 i i; n i 
ni,ki=even x z x ^ 



1 /• d^fc 



/ 7^^^''^'''' ^'' ^^ + ^^''^'''' ^'' ^))^r(^2, A;2; p) . (3.4) 



(p + n2qi + ^2^2)^ - mi 7 {2tt) 
Eqs. (|33), (^, (U) and (jS^) look very similar to the ones ^ for the 0(A^)-case. 



To rewrite Eqs. ( |3.1| ), ( p.2|) , (|3.3| ) and (|3.4| ) in a more convenient form, let us introduce 
the generating functions 

4«(.) = ™.Se".' + E^ ^,,^^,„^^!^^f._„^/ """*'""''« (3.5) 



*Here and below the poles should be understood according to the Feynman prescription rrr^ -^ rrr^ — iO. 
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7 ^"'^ 1 

and 

^^"'^^'^) = ^^+ ^ <(n,fc;p)J^e("-+'=-)-mX. (3-8) 

■^ 2 n,k=even 

Eqs. ( |3.1j ), ( p.2|) , ( p.3|) and ( |3.4| ) can then be rewritten, respectively, as 

||^-m?-t;(r)|$'^(r)=0, (3.9) 

||^-m^-^(r)|vl/"(r) = (3.10) 

||^-^?-t^(r)|e^-i-D;^V;P)=e^\ (3.11) 

||^-^2-^w} e-i-Z)5V;P)= e- (3.12) 



and 



where e is the energy component oi p, p = (e,p). 



^^1,2 = yp^ + "^1,2 (3-13) 

and 

t;(r) = A($^(r) + ^^^(r)) + >^ j ^S^^'^P) + ^f (^^P)) • (3-14) 

The next step in the transformation of Eqs. (U), (|3lOD , (|3Tl| ), ( ^l2D and (|33^) is 
to put D"'' and D" in the mixed representation — the coordinate in (imaginary) time 
and momentum for space. One defines 



D^'iu;,;T,r')=l^e<-^'^DfiT;p), (3.15) 



de 

2^ 

de 



Df{u,;r,r')=J^e^^-^'^Df{T-p). (3.16) 

In order to see that this quantity is indeed associated with the Fourier transform of 

amphtudes in energy, we notice that the insertion of Eq. (|3.7| ) on the r.h.s. of Eq. (|3.15|) 

yields 

r df: °° 

Df (^i;r,r') = llj- ^(^+nm,+km,)r-er' ijo, ^n; p) . (3.17) 

One recognizes now that e + nrrii + km2 is the energy component of the 4-momentum 
p + nqi + kq2 of the incoming particle with the 0(A'^i)-index a while e is that of p for b. 
In particular, the free propagator in the mixed representation is given by 

Dl\u,■,r,r')=^6^' / -^ , / . .^, = r^— e"-!— I (for A = 0) . (3.18) 
^ J 27r (e^ — a;i + zO) 2a;i 

11 



The equations ( p.ll| ), (|3.12|) and ( p.l4|) can be finally rewritten in the mixed represen- 
tation as follows: 

|^-o;^-^(r)|i-DjV;r,r') = -5(r-r'), (3.19) 

{^2-^"-^^)] ]^<(^' ^' ^') = -'^(^ - ^') (3.20) 



and 



v{t) = A(<l>2(r) + vl/2(r)) + -^ f / du uJu^ - miD'^{uj- r, r) + 

ZTT \Jmi 

+ f duuj^uj'^-mlDf{uj]T,T)] . (3.21) 



im2 

■\bh(, ,. ^ ^l\ „„A nl^Pi 



To understand the meaning of the summed amplitudes D^{uj; r, r') and D^ {oo; r, r') 

be 

X 



let us note that D^^{uj; t,t') has in the index space the structure 



(:bcc\ cb(:c 

^IG^(u;,r,r) + ^ 



D'^iu; r, r') = U''^ - ^ GJ(u;; r, r') + ^G^(a;; r, r') . (3.22) 



The amplitudes G^ and G"^ are associated, respectively, with the tensor and singlet 
0(iVi)-states of the two incoming particles. The averaged over the 0(A''i)-indices of two 
incoming particles amplitude is 

while the generating function for the symmetrized over all the n+2 0(A^i)-indices ampli- 
tudes is given by G^: 



each 



D^^(u;;r,r')^ = G^(c.;r,r'). (3.24) 



One sees that the averaged amplitude which enters Eq. ( |3.21|) coincides with G^ at large 

iVi. 

We can now rewrite Eqs. ( |3.19 ) and (|3.20|) at large A^ as the following equations for 
G^(a;; r, r') and G^{uj; r, r') (the latter is defined similarly to Eq. ( p.23| ): 

[^.-^"- <-')] G'x(^; ^' ^') = -^(^ - ^') ' (3-25) 

1 1^ - C.2 - t;(r) I Gj(o.; r, r') = -8{r - r') , (3.26) 

while v{t) is related to G'^{uj; r, r) and G^ (cu; r, r) by 

t;(r) = A($2(r) + ^\r)) + ^ f duj coJu^ - m? Gj(cu; r, r) + 

ZTT Jmi 

+ ^ / rfu; co^cu^ - ml GJ(cu; r, r) . (3.27) 

12 



Therefore, only the tensor amphtudes enter at large A^i,2- 

The transformation from Eqs. Q, Q, (jSj) and Q into Eqs. (gj), {^J§, 
(^), (^) and ( 13:271) is similar to the one of Refs. [0, ^, |15|, |0|, 0. These equations 
form the closed set which will be solved in the next section. 



3.2 The operator derivation 



The structure of Eqs. (|3.9| ), (|3.10|) , ( p.25| ), (|3.26|) and (|3.27| ) is easy to understand through 
the use of the functional technique |jTO| which relates r to the time variable t = xq oi the 
fields x"'{^) cind 0"(x) by 

T = it. (3.28) 

Using the definitions (|2.2|) and (|2.3|) , one rewrites ( ^3|) and ( |3.6|) as 

+ E 4i:T(^i---&n,A.../?fcix"(o)|o)ei'^---ei'"e2^^---e2^'=e("™i+'='"^)- 



1 



n,k>l 



n\k\ 



Y: ^TT(&l...&n,/3l...A|x"(0,-*r)|0)e?^---eH2^^---e2^ 



(3.29) 



and 



^°(r) 



m2^^e 

+ E 

n,fc>3 



« m2T 



n!A;! 



bn,(3i...(3k\r{omei'---ei"^^'---^^'e 



''1 . . . C^n^Pl . . . tA p {nmi+km2)7 



n,k>l 



n\k\ 



J2 {b,. . .6„,/3i . . ./?,|r(0, -^T)\0)e^' ■ ■■ei"^^2' ■ ■■^2 



(3.30) 



Similarly, using Eq. (|3.7|) and (13^^ ), one sees that v{t) which is defined by Eq. (|3.14|) is 
nothing but the sum of the matrix elements: 

v{t) = 
^ E ^7:1 (^1 • • • ^- A • • • /5.I : (X'(0) + 0^(0)) : |o) ^t^ ■ ■ ■ ^1-^2^^ ■ ■ ■ ^2^'= e (--1+^-^)- 



n,k=0 

00 



00 1 

^ E ^(&i...&n,A.../5.|:(x'(O,-^r)+02(o,-^r)):|o)et^---eH2^^---e2^n3.31) 



n,fc=0 

where : ... : stands for a normal odering. 

To simplify the notations, let us introduce the coherent state 

00 -| 



(3.32) 
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Then Eqs. (|3.29|) , (|3.30|) and ( p.31|) can be written as follows: 

$"(r) = (s|x"(0,-^r)|0; 



v[/-(r) 



'(0,-^r)|0 



and 



vir) = (E\: {x\0, -it) + 0^(0, -tr)) : |0) . 



(3.33) 
(3.34) 



(3.35) 



Eqs. (pD, ( glOD , ( g:25D , i ^M ) and {^^) can alternatively be deduced directly from 
the operator equations for the quantum fields x" and 0": 



and 



(-9^9^ - rnDr - A : (x' + 0')r : = 0. 



(3.36) 



(3.37) 



Let us take the matrix elements of Eqs. (|3.36|) and ( p.37|) between the states (S| and 
|0). We get 



+ 7T^-mn(S|x"(f,-^r)|0) 



and 



-A (H| : (x'(f , -tr) + 02(f , -2r))x"(f , -^r) : |o) = 



n2 Q2 ^ 

+ ^-ml\{E\r{S,-^r)\0) 



(3.38) 



[ 9r^ dap 

-A (S| : (x'(x, -tr) + 4>\x, -zr))0"(x, -zr) : |o) = 

Due to translational invariance one has 



(3.39) 



(3.40) 
(3.41) 



(H|x'^(x,-zr)|0) = $'^(r), 

(S|r(5^,-^r)|0) = vl/-(r) 
so that the partial derivative w.r.t. x vanishes. 

In the large A^ limit one can split the matrix elements of the operators : (x^ + 0^)x" : 
and : (x^ + 0^)0" : as follows 

(ai...a„,A.../3,|:(x' + 0')x":|o) = 
X (p(a„,+i) . . .p(a„),p'(/?fe,+i)) . . .p'(/3fc)|xlO) , (3.42) 
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where p, p' stand for permutations. This formula holds for n and k finite as N ^ oo and 
extends the standard factorization of 0(A^)-singlet operators at large A^ (see e.g. Ref. p9| ) 
to the case of multiparticle production. 



Using the definitions (|33^ ), ( [OSD , {^^ and ( p35| ), we rewrite Eq. (|^ in the 
form 

(e\ : (x'(x, -2r) + ^\x, -«r))x"(x, -«r) : |o) 
= (S| :(;^2(^^_^^)^^2^^^_^^)). |o)(S|x'^(f,-2r)|0) =t;(r)<l>"(r) (3.43) 

and analogously 

(S| : (x'(f, -it) + 02(f , -^r))0"(f , -«r) : |o) 
= (h| : {x\x, -it) + c^\x, -it)) : |o) (S|r (f , -zr)|0) = v{t)^^{t) (3.44) 

where translational invariance of the averages, which results in 

(h| : {x\x, -it) + 02(f , -tr)) : |o) = v{t) , (3.45) 

and Eqs. ( ^40]) , ( Fill) have been used. Substituting ( ^4^ ) and ( ^441 ) into Eqs. (|3]3|) 
and (13:^ ), we obtain Eqs. (|J) and ( ^lOD . 

To infer the remaining Eqs. ( |3.25|) , (|3.26| ) and (|3.27|) , let us represent D'^{uji]t,t') 
and D'^^{uj2] t, t') given by ( |3.15|) and (|3.16| ) as 

d?x 



Df{oo,;T,T') 



(27r) 



and 



Df{u;,;T,T') = J 



<P 



X 



(2vr) = 



^e'P^E\Tx''{x,-tr)x'{0,-tT')\0 



-S'^^\p)<^%t)<^\t') 

S|T0"(f,-ir)0^(O,-zr')|O 
-5(=^)(p)^"(r)*^(r'; 



(3.46) 



ipa; 



(3.47) 



Here the subtractions cancel disconnected parts which are not included in the definition of 
Z^"^ and D2 and T stands for the T-product. Notice that Eq. (|3.27|) can immediately be 
rewritten in the form ( p.35 ) using ( 3.46|) and ( 3.47|) . The terms involving the propagators 
D""^ and D°"^ in the equation for v{t) are associated with the connected part of the 



operator : x^ + 0^ : 



The operator equations ( p.36| ) and (|3.371 ) result in the following equations for the 
T-products: 



d^ 



+ 



d^ 



^Q^. ' Q^ m?|(H|Tx^(f,-.r)x^(0,-^r')|0) 

-A (s|T(x'(x, -it) + 02(f, -tT))x\x, -^r)/(0, -^r')|o) 

= -5^''6{t-t')6^''\x) 



(3.48) 
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and 



92 



92 



-A (s|T(x'(x, -it) + 02(f, -zr))0"(f, -zr)0^(O, -^r')|0 



(3.49) 



where the delta functions on the r.h.s.'s emerge in a standard way when differentiating 
the T-product and using canonical commutation relations. 

Eqs. ( p.48| ) and ( p.49| ) are the lowest ones in the chain of the Schwinger-Dyson equa- 
tions which can be extracted from the following identities 



and 










(50° (x) 



■F[xA] 



i{ ^ 




(3.50) 



(3.51) 



where S is the action of the model and F[xi 0] is an arbitrary functional of x ^"^^ (p- 
In particular, Eqs. ( p.38| ) and ( |3.39| ) are associated with F = 1 while Eqs. (|3.48| ) and 
( p.49| ) correspond to F = x^(0, —ir') and F = 0'^(O, —ir'), respectively. In the language 
of path integral Eqs. ( |3.50| ) and ( |3.51D result from the invariance of the measure under an 
arbitrary variation of x°'(x) or 0°(x). 

Using the large- A^ factorization 

[e\T{x\x, -zr) + 02(f, -^r))x'^(x, -^r)/(0, ~it')\o) 



= {e\ :(x\x,~tT) + <p\x,-tT)): \o) {E\Tx%x,-tT)x\0,-tT')\o) 

and analogously 

(s|T(x'(x, -zr) + 02(f, -zr))0"(f, -zr)0^(O, -2r')|0 
I ■.(x\x,-ZT) + ^\x,-tT)): |o)(s|T0"(x,-zr)0^(O,-2r')|o) , 



(3.52) 



(3.53) 



as well as Eqs. ( p.45|) , ( p.9|) , ( p.lO|) and upon integrating over d^x, we obtain Eqs. ( p.25|) , 
(g) from Eqs. dOHQ , (|^. 

Some comments concerning the operator derivation of Eqs. (|3.9|), (|3.10|) , ( p.25| ), ( |3.26|) 
and ( |3.27| ) are in order. The factorization was crucial to drop out the dependence on the 
spatial coordinates in the interaction terms and, therefore, for the reduction at large N 
to the quantum mechanical problem which is given by the closed set of Eqs. ( |3.9| ), ( |3.10| ), 
( p.25| ), ( p. 261) and ( p. 27]) Its appearance is due to the special kinematics of the produced 
particles whose wave functions do not depend on the spatial coordinates. The factorization 
formula (|3.42| ) holds for n and k fixed as N —>■ 00 (notice that only ni 2 < n and ki^2 < k 
enter Eq. ( |3.42| )). The corrections to the large- iV amplitudes are known [0, [1^ to be 
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controlled at large n by the parameter n^X ~ n^ /N (also by k'^ /N in the given case) 
and are due to rescatterings of the produced particles. Hence, equations (|3.9|), ( p.lO| ), 
( p.25| ), ( |3.26| ) and ( |3.27| ) allow for the calculation of the multiparticle amplitudes only for 



4 The exact solution 

To solve the set of equations ( |3.9| ), ( |3.1CI| ), ( |3.25| ), ( |3.26| ) and ( |3.27| ), let us first look at 
Eqs. ( |3.25D and ( ^.2(j| ) for a given v{t). These equations determine the Green functions of 



the Schrodinger operator with the potential v{t) while r plays the role of a 1-dimensional 
coordinate. Equivalently G'^{uj] r, r') is the matrix element of the resolvent 

G^^-.ry) = {r\ ^,1 ,^ \r') (4.1) 

where D stands for d/dr for brevity. One should then take the diagonal matrix element 
of the resolvent, C^iu] r, r), to substitute it into Eq. (|3.27|) and determine f (r) versus ^^ 
and \l/^. 

The general solution of this problem for arbitrary v is given by the Gelfand-Dikii 



formula 30 



G>;r,r) = i?.H^f:^ (4.2) 

where the differential polynomials Ri[v] are determined recurrently by 

Ri[v] = j,(^Id'-v-D-\D~^' .^- (4.3) 

and the inverse operator is 

D-h{T) = r dx v{x) . (4.4) 



Eq. ( [4. 3D stems from the fact that Ruj[v] obeys the third order linear differential equation 

i (1^3 -Dv- vD) R^[v] = u^DR^v] . (4.5) 

The polynomials Ri[v] depend on v and its derivatives f*^*-* = {D^v). The first few poly- 
nomials are 

RM = \. Ri[v] = -\. R2[v] = ^{"iv' - v") , 

RM = - — (10t;3 - low - h{v'f + t;(^)) , . . . (4.6) 

64 

while for r-independent v{t) = vo one has 

1 

Ru,[vo] = o / 9 I ( Wo = const.) (4.7) 
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which agrees with Eq. (|3.18|) at r = r'. 

Following Ref. ^6\ we renormalize the above equations. Let us introduce Vq as 
A 



-^0 



Ni f doj^Juj'^ - mlj^ + N2 f du 



UJ" 



m. 



2R I 5 



4:71 \ J mm ''m2R 

where the renormalized masses mi,2R are defined by the following equations 



m 



IR 



ml + vo, ml^ 



ml + vq. 



We also introduce the renormalized coupling constant defined as 



1 

A^ 



+ 



duo 



^ ^-^Ir N2 



A 87r^ JmiR ijJ 

The renormalized potential f_R('r) reads 

VR{,r) = v{t) - Vq 



duo 



u" 



m 



2R 



TX Jm.2R 



UO^ 



(4.8) 



(4.9) 



(4.10) 



ZTT JmiR \ LO 

+iV2 / duo uo^uo^ - ml^ [r^Ivr] - ^^ 

Jm2R \ UO 



2uo 



2uo 



Now everything is expressed in terms of mi^2R and Vr{t): 



(4.11) 

(4.12) 
(4.13) 

^ = ^Jp^ + ml2R. (4.14) 

The integrals over uo on the r.h.s. become convergent after the renormalizations. The 
meaning of renormalization is that one chooses the bare quantities, ml 2 and A, to be 
dependent on the cut-off according to Eqs. ( [4.9|) and ([4.1 0|) and render the renormalized 
ones, m|j and Aj^:, cut-off-independent. 

We can easily find an exact solution to these equations. Let us denote 

L = 1 (^ID^ -Dv- vD) . 
We shall look for R^ [vr\ in the following form 

R^[vr] = a{uj) + x(cu)$2 + y{uj)^^ 



where a, x and y are constants which may depend on uo. Using eq. 
that $^ and \E'^ are the eigenvectors^ of the operator D~^L : 



(4.15) 

(4.16) 
and the fact [BHl 



1 (iZ}2 _ ^^ _ D-'vnD) $2 



"i2 „. 7~i — 1„. 7~i\ /^2 1 n^ 



^ ' iD^ - Dvr - vrD) $2 = D-^Dml^^'^ = m?^$2 



(4.17) 

'''This is derived using solely Eqs. (4.12) and (4.13) for an arbitrary u_r(t) and the fact that $^(t) and 



\I' (r) vanish for r — > —00. 
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I [ID^ -VR- D-\rD) ^2 = iZ)-l (1^3 _ ^^^ _ ^^^^ ^2 ^ 2?-l/5^2^$2 ^ ^2^^2 ^ 

(4.18) 
one gets 



LRu^[vr] = --Dvr + xmlj^D^^ + ymlj^D"^^ 
We remove D from both sides of the above equation and we get 



(4.19) 



2x 



2 , .2\^2 I ^y /^2 , ,2\,T,2 



Vr = c+ —{ml^ - uj')<l>^ + —{m^R - ^')^ 
where c is a constant of integration. Thus we have 

vr = c + a$^ + /?^^ 
where the constants 



a 



2x{u? - m\j^ _2y{uPj-n4^ 



(4.20) 



(4.21) 



(4.22) 



do not depend on u since vr does not. 

We now substitute both expressions for R^ and vr into the equation (|4.11| ) which 
determines vr. One can easily see that c = 0, a = 1/(20;) (notice that this value of a 
agrees with the perturbation expansion) and 



A 



a 



P 



R 



a 



1 + ^A' 



X 



/3 

A^r ' y A,,(,a __2 V 

87r2-'2 



87r2- 

\r 



Auj{uj^-mlj^)' 



where 



2 
IR 



h = Ni dujyuj^ — m 

+N2 / dcuJcu^ — m. 



uj^ — m 



IR 



uo^ 



+ 



2 

2R 



■'n-2R 



ijj'- — m 



IR 



LU^ 



(4.23) 
(4.24) 



(4.25) 



IR 



l2 = Ni f duj^uj"^ - ml 

JrniR 

+N2 / dusjcj'^ — ml 

Jm2R 

After some calculations one gets 



uj'^ — m,2R uj"^ 



2R 



u"^ — m,2R uj"^ 



I^= Ni + N2 f duj^Juj'^ - m 

Jm2R 



2 
2R 



uj'^ — m\^ u;2 



(4.26) 



(4.27) 
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and 

1 1 



JmiR 



(4.28) 



For nil = "^2 these formulas recover the resuhs |26| for the 0(A'^i+A^2)-case. The fact 
that the solution is unique can be proven quite similarly to Ref. [|^. 

The diagonal resolvent reads 

R4vr] = 7r- TT-2 —^ - TT^ ^- 4-29 

The appearance of the poles in the resolvent has two faces. 

1. From one point of view it is not directly related to the propagation of particles in 
the loops. For example the pole in the resolvent appears for the symmetric 0{N) case even 
at the tree level. The pole is actually due to an existence of a strong classical external field 
which induces a potential v{t) and corresponds to a bound state of a quantum mechanical 
particle. Thus the pole in the resolvent is actually a manifestation of a bound state in 
the spectrum of the corresponding Schrodinger operator. This is of course an obvious 
fact but it explains the appearance of an extra pole when we switch on an additional 
particle. Emissions of an additional particle induce a new term in the effective potential 
which is not suppressed at large N since we assume that the wave function of an emitted 
particle is of order a/ZV. This new potential is responsible for a new bound state of the 
quantum mechanical particle. This also explains why the residue of the additional pole 
in the resolvent is not suppressed at large N. 

A natural question would be as to why the poles should be exactly at the values of the 
masses of the particles (a priori these poles could be anywhere). The answer is simple: 
this is because we know that there exist solutions [^ to the Schrodinger equation for 
uj = mi^2 which decrease at both infinities. Thus we know that indeed they are the bound 
states with the energies mi^2- 

2. On the other hand the diagonal resolvent itself corresponds to a loop by definition 
since it corresponds to propagator integrated over the energy {i.e. to the one with coin- 
ciding ends in the r-space). For simplicity let us consider the model of the 0(A'^i)-particle 
with the mass mi and an additional singlet one with the mass m2, i.e. N2 = 1. In per- 
turbation theory the simplest correction to the propagator of the 0(A'^i)-particle is given 
by a diagram with a 4- vertex which produces 2 outgoing particles with the mass m2- To 
extract the diagonal resolvent we should make a loop of this propagator which is done in 
the momentum space by integrating over the energy (de) of the virtual particle. Thus we 
easily see that the correction to the resolvent which is due to an emission of two particles 
is proportional to 

.00 2^ (c^2 _ ^2 _ ,o)(cu2 - (e - 2m2)2 - zO) ' ^^'^^^ 

where u'^ = |pp + ml, and p is a spatial virtual momentum. This integral is convergent 
and can be computed by closing the contour down around e = u — iO and e = uj — 2m2 — i0. 
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It is easy to see that a pole appears at c<j = m2: 



I 



^L'^ oc — r. (4.31) 

Am2Uj{uo - m2 - iQ) 

This is exactly what is implied by (|4.29|) . 

We thus have checked that there should be two poles in the resolvent in the model 
of two particles with different masses. These poles sit exactly at the masses of those 
particles. 

We now can find the expressions for $" and \1'". We have the following equations 

[d^ - m\^ - (a<l>2 + /3^2) I $a _ 0, (4.32) 

These equations reduce to the classical ones (with the renormalized coupling constant A/?) 
if we renormalize the fields $° and \1/" as follows 

^-= fK^-=ll + h^l\ $-, ^,-= lKm-=(i + h^i\ ^". (4.34) 

V Ar V Stt^ y \ Xr \ stt^ ) 

The solution to these equations is given by the expressions which are a straightforward 



generalization of those found in Ref. P3| : 



^^ =~zl{l- 2XR^~zi\ (l - ^Sf - ^~zi + Xl^^-zf-zl] ' , (4.35) 

^^ = t f 1 + 2A«4-^il f 1 - ^^i' - ^^2 + >^\-^^~zl~^^ ' , (4.36) 



where 



and 



miR - m2R 
miR + m2R '' 



/ \ \ -1/2 



(4.37) 



~A = ^2 ^T^e'"-^ = es^e'"-^ 1 + -4/2 . (4.39) 

V Ar \ stt y 

Notice that the constants a and /3 are in general complex. By assuming for definiteness 
that m|^ > mlji one can easily see that the constant a is real while the imaginary part 
of j3 is given by the following equation 

I„4 = -^ . \/§E^. (4.40) 

(3 8n 2m2R ^ ' 
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Figure 5: The one- loop diagram for the process x + X ~^ + </• which possesses an 
imaginary part for m2R > miR. given by Eq. ( 4.40| ). 



The appearance of this imaginary part is in contrast to the 0{N) model with unbroken 
symmetry where the amphtudes are real pGf. In the present model the imaginary part 
appears in a or /3 because the masses of 0{Ni) and 0(A^2) particles are not equal. Indeed 
in a particular example ^2/? > ''t^ir ^in imaginary part appears in the amplitude of 
scattering of two x particles due to a possibility of an inelastic process X + X ~^ + 0- 
This imaginary part is associated with the diagram of Fig. |^. A direct calculation of the 
imaginary part of the diagram of Fig. ^ gives the same expression ( |4.40| ). 



In the case of A^2 = 1 which is discussed above, we drop out the contributions pro- 
portional to Z)?" in the recurrence equations. Therefore this case can be easily recovered 
from Eqs. (|03| ) and (g]^ by taking formally A^2 = in Eqs. (g]^ and (g]^ for 



Ji and I2 which enter the expressions ( [4.23| ) and ( |4.24| ) for the constants a and (3. The 
correct imaginary part immediately appears for 1712 > mi after integrating (|4.31| ) over u. 



This imaginary part corresponds to on- mass-shell virtual 0{Ni) particles emitting heavier 
singlet particles (the wave functions for incoming particles are factorized out from this 
diagram and therefore it can be thought as a contribution for the forward scattering for 
singlet particles). 

A generalization of the equations for the generating functions $, \E' to the case of 
0{Ni) X ... X 0{Ns) scalar field is straightforward (s is a positive integer). In such a 
general case which is described by the Lagrangian 

^ = lti9,<t>^)i9,<t>^)-lt^"<l^"-jit<l^^ ' (4-41) 

the equations for the generating functions for the amplitudes of the multiparticle produc- 
tion at the threshold by a single particle are given by equations similar to the classical 
ones while the effect of loops in the large N is reduced to a renormalization of the masses 
and the coupling constant and to a (complex) renormalization of the wave functions ^i of 
the external particles. 

The potential of the associated quantum mechanical problem is 

VR = Y.a,<^l (4.42) 
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where $f are the generating functions for the amphtudes of multiparticle production by 
a single particle 0°, and the parameters a^ are defined as follows 






UJ^ — m^^ uj 



(4.43) 



Here rriifi stand for a renormalized mass of the i-particle. The diagonal resolvent reads 

This resolvent corresponds to a nullification of all amplitudes except for those associated 
with the production of two identical particles. 



5 A model with spontaneously broken symmetry 

5.1 The Schwinger— Dyson equations 

Let us now consider the model of an 0{Ni) scalar field and a singlet scalar field with the 
wrong sign of the mass square. The Lagrangian reads 

£ = i{d,x'){d,x') + U9Md,<i^) - X(xY) - 1^2(0') - -A<P<P + x\')\ (5.1) 

where m\ < and ml > —Im^l can be either positive or negative Q. The interaction term 
has an 0(iV)-symmetry with A^ = A'^i + l which is explicitly broken for ml 7^ m"^ by the 
mass term. The Lagrangian (|5.1|) is invariant for ml 7^ m^ under the reflection: -^ —0. 

At the classical level one can see that the reflection symmetry is spontaneously broken 
and the field acquires a vacuum expectation value 



F^2 



while the masses of the physical particles are 



(5.2) 



m.^ = J mil + I'^il -I ""^(A = v2|^7i2| (tree level) . (5-3) 

If the 0(A^)-symmetry were not broken explicitly by the mass term {i.e. if mi = 7712), one 
would get A^ — 1 Goldstone bosons with the mass m,^ = 0. 

The classical equations of motion read 

{5«-|^2l + A(xV + 0')}0 = O. (5.4) 



'Notice that in this case the vacuum state corresponds to an unbroken 0{Ni) symmetry. 
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At the classical level one can shift the field = 0' + 0o- Then the classical equations 
of motion take the form 

{dl + {ml + \ml\) + 2A0O0' + Kx^x'' + <P")}x' = , 
dl<P' + 2\ml\<p' + 3A0O0" + A^oX'^x" + A(xV + <P"W = 0- (5.5) 

One would expect that Eqs. (|0| ) get simply generahzed to the quantum level through the 
substitution of a normal ordering into (|5.5|) . However this is not quite correct because the 
vacuum expectation value 0o can be renormalized due to the (p'x^ interaction. Therefore 
it is convenient to consider the equations of motion for the unshifted field 0. 

One can extract the equations for the amplitudes to produce physical particles by 
using the same method with the case of no spontaneous symmetry breaking. We shall 
assume that the field contains a constant part which corresponds to a correct (at the 
quantum level) vacuum expectation value of 0. We then take the matrix elements which 
are associated with the production of physical particles at the threshold when the wave 
functions of the produced particles are translationally invariant similar to the procedure 
followed for the case without spontaneous breaking of the 0{N) symmetry. 

Indeed we have the following operator equations 

id,d, + rnDx"" + X:ix' + 0')X'^ := (5.6) 

and 

{d,d, - \ml\)<f) + X:{x' + </>')</> := 0, (5.7) 

where : : stands for a normal odering. Let us take the matrix elements of these equations 
between the state (ai . . . an, 0'^| and the vacuum where 0' stands for a physical particle 
that corresponds to the field 0. 

One can see that 

^'(^)= E -^(«i---«n,0'1x'^(O,-^r)|o)eJ^---eiH2' (5.8) 



,^-1 n\k\ 



and 



^M= E -^Aai...an,<P"\<P{0,-zT)\0)ei'---ei"e2. (5.9) 

n,k>l "'•'^• 

are nothing but the generating functions for the multiparticle production from a single x 
or 0' particle, respectively. By introducing again the coherent state 

(S'l = E -jn^l' ■ ■ ■ ^1-^2 («1 • • • «n, 0'1 , (5.10) 

n,k>l 



we rewrite Eqs. (|5.8| ) and (|5.9|) as 

$«(r) = (H'|x'^(0,-^r)|0), (5.11) 
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vl/(r) = (H'|0(O,-2r)|o). (5.12) 

In the large A^ limit one can split the operators : (x^ + 0^)x" • cind : (x^ + 0^)0 : as 
follows 

(oi...o,„.*"|:(x' + .*V:|0) = 

= 1: E -^;^(p(a0...pK).*'"|:(X^ + *=):|0)x 

x(p(a„,+i)...p(aO,</)"1xl0) (5.13) 

and 

(ai...a„,0"=|:(x' + 0')0:|o) = 

x(p(a„,+i)...p(a„),</>'^^|</>|0) (5.14) 

where p stand for permutations. Using the definition (|5.1CI| ), Eqs. ( ^.13| ) and (|5.14|) can 
be rewritten in the simple form 

(S'l : ix' + 0')X" : |0) = (s'| : (x' + 0') : |o) (S'|xl0) (5.15) 



and 



(S'l : ix' + 0^)0: |0) = {e'\ : {x' + <P') : |o) (SV|0) . (5.16) 



The Schwinger-Dyson equations can now be derived quite similarly to Subsect. 3.2. 
It is crucial again that one drops out the dependence on 3D spatial coordinates in the 
interaction terms. This is correct due to the factorization of these terms in the large- A^ 
limit and the special kinematics of the produced particles. We thus get the following 
equations 

{D''-ml-v{T)}x' = 0, (5.17) 

{d^ - v{r)} (p = 0, (5.18) 

where 

v{r) = A($^ + ^'- <Pl) + V (0i^x (^;P)- (5-19) 

Here we have taken into account that |m2| = A0q and used the physical mass m^ of the 
field X which appears when we take into account the non-vanishing vacuum expectation 
value of the field (p. At the classical level m^ is given by Eq. ( |5.3| ). 

The propagator V^ obeys Eq. ( p.ll| ) as for the case without spontaneous breaking. 
Introducing D'^{uj; r, r') (the mixed representation) as in Eq. ( p.l5| ) obeying 

{Z^2 _ ^2 _ ^(^)| l-D'^ico; T, r') = ~5{t - r'), (5.20) 
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where iJ^ = rn^ + |p|^ (p stands for the spatial components of momentum), we can rewrite 
the potential v{t) given by ( p.l9| ) as follows 



^;(r) = A(<l>2 + vl/2 _ 02) ^ AiV r ^^^ /^2 _ ^2G'^(^. ^^ ^)^ (5.21) 

ZTX Jrrix 

where G"^ stands for the tensor part of D'^ . Actually these equations can be easily obtained 
through the use of generalized recurrence equations for the amplitudes to produce particles 
X"" and (f)' . To do that we have to use the Lagrangian for the fields x ^i-iid (j)' = (p — cj)^ which 
are physical at the classical level. The only difference with the case of no spontaneous 
breaking is a presence of the vertices x^0' and 0'^. 

We now renormalize the masses and the coupling constant. Let us introduce 



The renormalized masses are then defined as follows: 

'^li? = m\ + vo, mlj^ = m\ + vq, (5.23) 

while for the renormalized coupling constant we have 



- = y + — / du^ ^. (5.24) 

We see that the physical mass of the field x is not renormalized, m^ = m\j^ — m\j^ = rn^. 
The renormalized potential is defined as 

vr{t)=v{t)~vq= (5.25) 



A^($2 + vl/2 _ 02) ^ _|^^ I ^^ ^^^2 _ ^2^ [R^[vr] 
zvr Jm^ B \ 



Here (f)Q = \itl2r\/^/Xr is the renormalized vacuum expectation value of the field 0. 
Now everything is expressed in terms of rrii,2_R and Vr{t): 

[D^-ml-vn{T)]^\T) = Q, (5.26) 



uo= ^p^ + ml. (5.27) 

The integral over uj on the r.h.s. of Eq. ( |5.25| ) becomes convergent after the renormaliza- 
tions. 
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5.2 The exact solution 



We shall look for Ru)[vb\ in the following form 

Ru.[vr] = aiu) + xico)^^ + y(^)^^ (5.28) 

where a, x and y are constants which may depend on u. Similarly with the previous case 



of no spontaneous breaking we shall use Eq. (O) and the fact that 



L$2 = 1 Ud^ _ Dvr - vrD) $2 = mjD$2 



L^^ 



2 V2 



./}3 _ ^^^ -VrD]^ 



We obtain 



LRu,{vr] = --Dvr + xmlD^'' 



(5.29) 
(5.30) 

(5.31) 



= uj^D{a + x(^^ + y^^)- 
By removing D from both sides of the above equation we get 

a ^ a 

where c is a constant of integration. Thus we have 

vr = c + a$^ + /3^^ 

where the constants 



a 



2x{u^ - m^) 2yu' 

5 P 



(5.32) 

(5.33) 
(5.34) 



a a 

do not depend on u since Vr does not. 

We now substitute the expressions for R^ and vr into the equation that defines vr. 
We have 

c + a^^ + (3^^ = XR{^^ + ^^-(f)l)+ (5.35) 



+ 



27r2 



dujuj^loj'^ — m^ 



aa 



a — 



2(u;2 — m?,) 



$2_ 



[3a 

2J^ 



^2 + 



Auo^ 



[c + a$2 + /3^2N 



From the above equation we get 

uXrN 



a = \r + 



27r2 



dujijjJijj^ — m 



Au^ 2(cj2-m2^_i 






dujujjuj'^ — m 



1 



4cj3 2cj2 



One can also extract an equation for the constant part of Eq. ( ^.35| ) 






27r2 L dujuj^uj^-ml 



j3a ^ 1 2\ 1 



2a; 

(5.36) 
(5.37) 

(5.38) 
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a/302 

2^2 


1 

2w' 


UJ 





and hence 

We substitute this expression into Eq. ( |5.38D we get 



We get three equations for four unknown parameters a, /5, a and c. 

To choose the correct solution to Eqs. ( |5.36| ), ( p.37| ) and (|5.38|) we have to use the 
following physical condition. When the resolvent is expressed in terms of the physical 
(shifted) fields ^°' and \1/' = \1/ — 0o the constant part of the resolvent is fixed in pertur- 
bation theory to be equal to l/(2u;). This means that 



(5.39) 
(5.40) 

(5.41) 

(5.42) 
(5.43) 

(5.44) 
(5.45) 



For the parameters a and (3 we have 



A; 



a 



/3 



1 + 



87r2 



-Jl 



X 



R 



XrN J ' 



Stt 



— •J2 



where 



Jl = / dujJuj^ — m 



U 



J,. 



dujJuj'^ — m 



(cj2-/502/2)(cj2_^2) ^2 

1 1 



_uj^-(3<Pl/2 uj\ 

It is easy to see that the parameters a and (3 are real in contrast to the case without the 
spontaneous breaking of the 0(A^)-symmetry. 

For the resolvent we get the following expression 

-I acj$2 /5^2 



Ru^[Vr] 



(5.46) 



2u{u^ - l3(t)l/2) 4(cu2 - m2/4)(cj2 _ ^2) 4^(^2 _ ^2/4) ' 

where we identified the physical mass of the field \1/' with 

ml = 2(3ct>l (5.47) 

One can see that the resolvent has two poles at u; = m^ and u = m^/2. It can be easily 
checked in perturbation theory with respect to the external lines at the tree level for the 
resolvent that exactly these two poles appear to the lowest orders. Indeed the pole at 
uj = m^ corresponds to a diagram with two external legs of the x fi^ld while the pole 
at u; = m^j2 comes from the diagram with a single external leg. It is worth noticing 
that by looking at Eq. ( |5.46|) one might think that there is a pole at a; = m^/2 in the 
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amplitude to produce two x particles which does not of course appear in perturbation 
theory. It is easy to check however that such a pole is cancelled if we take into account 
the explicit solution for $'^ and ^ which is given below. One can also see that there are 
no poles at C(j = m^ in the resolvent. This signals the fact that the amplitudes for the 
production of two external particles are cancelled in contrast to the case without the 
spontaneous breaking. 

Thus we get the following equations for the generating functions 



[d"^ - (a$2 + /5^2 - /30^)} ^ = 0. 



(5.48) 
(5.49) 



These equations reduce to the classical ones (with the renormalized coupling constant A/j) 
if we renormalize the fields $" and \1/ as follows 



^'' = ,1 ^^'^ = ^''(1 + 4^^i 



-1/2 






^1/2 



57r" 



The solution to the above equations which obeys the condition that at A/? — i> 



(5.50) 



(5.51) 



is similar to that given in Ref. p3 



^Oi? 



1 



1 + 



Z2 



Z2 



+ 



8Ar 2 , 4A^^ 2m,, - m^ .\ 

r^i H ] 77; : T^^2^i X 



^Oi? 



Aml-ml^ (poR {2m^ + m^)^ \ 



(5.52) 



8Ar ^2 , 4Ak 2my, -m^ ^ ^^ 

-z^ H — ] — t:^ : r7^2^i 



-1 



20oi? 4m2-m^ ^ 0oij (2m;^ + m^)^' 



$'' = zM 1 



Here 



and 



2m^ — m^ ^2 
2m^ + m0 20oK, 



1 - 



Z2 



\i? 



^2 , 4Ar 2m^ - m^ ^ ^2 

2:2^1 



^O/J 



-^r + 



4m2 - ml (pon {2m^ + m^f 



%R 



\ Xr V 8^ / 



-1/2 



5i = m^Cie' 



Z2 = m<^^2e" 



A^ 

87r2 

AfiiV 
87r2 



-Ji 



-J2 



-1/2 



-1/2 



Thus we finally get 
^ = 0( 



1 + 



^2 



+ 



iX 



R 



?0R 



2 4Afi Jm^^m^ 2\ ^ 
^1 + 1 — To ; \^^22:i X 



Ami ~ ^l ^OR (2"^x + ^<p)^ 
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(5.53) 
(5.54) 

(5.55) 
(5.56) 

(5.57) 




Z-i H ; -— —Z2Z1 



^ \ 2m^ + m^ 2(j)Qji J \ 20oh 4m2 - m^ ^ 0oi? (2m^ + m^)3 ^ ^ 
We observe that this solution is real. 



(5.5^ 



6 Conclusions 

Let us briefly summarize our results. We presented an exact solution in the large- A^ limit 
for the amplitudes of multiparticle production at threshold in the 0^ theory with the 
0(A^)-symmetry which is softly broken to O(iVi) x . . . x 0{Ns) {N=Ni+ . . . +Ns) by the 
mass term. We found that the effect of loops in this limit reduces to a renormalization of 
coupling constants and masses. Moreover as the masses rrii of the 0(iVj)-multiplets are 
nonequal, the amplitudes of multiparticle production become complex. The explicit form 
of the solution demonstrates a nullification of the multiparticle amplitudes with all legs 
on mass-shell. The only non-vanishing amplitudes correspond to the processes 2— i'2 when 
the incoming (outgoing) particles are identical. 

We also analyzed the model of the 0{N) + singlet scalar particles in the case when the 
reflection symmetry is spontaneously broken. In the limit of large A^ the effect of loops 
of the 0{N) field reduces to a renormalization of the coupling constants and masses of 
the physical particles while no imaginary part of amplitudes appears in this case. Using 
the exact solution for this problem, we demonstrate a nullification of amplitudes except 
for those with 2 incoming 0{N) particles and 2 outgoing 0{N) particles or 1 outgoing 
physical singlet particle. 



Some problems remain unsolved. As it has been demonstrated in Ref. [^, the nulli- 
fication of the tree level amplitudes for the model of two scalar fields with further restric- 
tions on kinematics, when one particle, say m2, can decay at rest into the others with 
the mass mi, is due to a symmetry of the classical Lagrangian where the dependence of 
the fields on the spatial 3D coordinate is discarded. It is straightforward to check that 
such a symmetry exists at the tree level in the 0{Ni) x 0{N2) model as well. Indeed 
the Lagrangian (^TT|) is invariant (up to terms which are total derivatives) for the fields 
independent of the spatial coordinates under the following transformations 

50° = e[-Ax"(x>" - rx") + \{rnl - m^)0«], 
where stands for d(f)/dt and e is a constant parameter of the variation. Following 



Ref. 1 24 1 one can see that the no n- vanishing amplitudes come from the resonance terms 
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which appear in the iteration procedure for solving the classical equations of motion with 
the boundary conditions at A -^ 0: 

X^Cie-'""'', 0^^^"e-^-2t_ (6.2) 

The explicit form of our solution shows that such a symmetry survives at the multi-loop 
level (at large N) since the only effect of multi-loop corrections reduces to a renormaliza- 
tion of the coupling constants and masses. Obviously we come to the same conclusion in 
the case of a spontaneously broken reflection symmetry in the model of the 0{N) + singlet 
scalar particles. An interesting question about this symmetry is whether it can be ex- 
tended to arbitrary mi and 1712 when kinematics requires nonvanishing spatial momenta 
of incoming particles. 

As it is already mentioned in Subsect. 3.2, our results for the multiparticle amplitudes 
are applicable when the number of produced particles is much less than a/ZV. They nicely 
show how the loop effects renormalize tree amplitudes but do not take into account dia- 
grams with rescatterings in the final state which are expected to restore unitarity. There 
exist at the moment two examples of such calculations: 1) a semiclassical calculation [^ 
for the A^ = 1 case with a spontaneously broken reflection symmetry, 2) models with any 
A^ in 2 + 1 dimensions which are solved [^] by summing infrared logarithms with the 



aid of renormahzation group. It would be interesting to develop a technique based on 
the Schwinger-Dyson equations which would be applicable for the case when the num- 
ber of produced particles is of order y/N. We hope that the operator formulation of the 
Schwinger-Dyson equations for the multiparticle amplitudes which is proposed in this 
paper may be useful for this purpose. 
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